Calculation of high- dimensional probability density functions of stochastically excited nonlinear mechanical systems

Technical systems are subjected to a variety of excitations that cannot generally be described in deterministic ways. External disturbances like wind gusts or road roughness as well as uncertainties in system parameters can be described by random variables, with statistical parameters identified through measurements, for instance. For general systems the statistical characteristics such as the probability density function (pdf) may be difficult to calculate. In addition to numerical simulation methods (Monte Carlo Simulations) there are differential equations for the pdf that can be solved to obtain such characteristics, most prominently the Fokker-Planck equation (FPE). W. Martens MMD TU Berlin Tel.: +49 030 314 22876 Fax: +49 030 314 21173 E-mail: [email protected] U. von Wagner MMD TU Berlin Tel.: +49 030 314 21169 Fax: +49 030 314 21173 E-mail: [email protected] V. Mehrmann Inst. f. Mathematik TU Berlin Tel.: +49 030 314 25736 Fax: +49 030 314 79706 E-mail: [email protected] A variety of different approaches for solving FPEs for nonlinear systems have been investigated in the last decades. Most of these are limited to considerably low dimensions to avoid high numerical costs due to the ”curse of dimension”. Problems of higher dimension, such as d = 6, have been solved only rarely. In this paper we present results for stationary pdfs of nonlinear mechanical systems with dimensions up to d = 10 using a Galerkinmethod which expands approximative solutions (weighting functions) into orthogonal polynomials.